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Abstract 

We consider the two-nucleon weak interaction with a pionless effective field theory. Dibaryon 
fields are introduced to facilitate calculations and ensure precision in the initial and final state prop- 
agators. Weak interactions are accounted for with the parity-violating dibaryon- nucleon-nucleon 
vertices, which contain unknown weak dibaryon-nucleon-nucleon coupling constants. We apply the 
model to the calculation of a parity-violating observable in the neutron-proton capture at thresh- 
old. Result is obtained up to the linear order in the unknown dibaryon-nucleon-nucleon coupling 
constants. We compare our result to the one obtained from a hybrid calculation, and discuss the 
extension to weak interactions in the few-body systems. 
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I. INTRODUCTION 



Weak nucleon-nucleon (NN) interaction has recently been formulated in the framework 
of effective field theory (EFT) [1.]. Parity- violating (PV) weak A A potentials have been 
obtained up to next-to-next-to leading order (NNLO) in the pionful theory. The weak 
potentials obtained from EFT have been subsequently applied to the calculation of PV 
observables in the two nucleon systems , and the results demonstrate the adequacy 

of perturbative scheme of the EFT for the description of weak A A interaction. 

About 40 years ago, Danilov suggested the parametrization of the parity mixing due to 
the weak A A interaction in terms of five PV low energy constants [3, Is L and the idea was 
applied to the calculation of PV observables in the few- nucleon systems 3, sj]. In this work, 
we will consider PV observables in the neutron-proton capture at threshold, where energy 
scale or momentum transfer is much small compared to the pion mass. At this small scale 
pion can be treated as a heavy degree of freedom. If pions are treated as heavy degrees, 
we can integrate out the pion from the theory, and then we obtain a pionless theory where 
interactions are described in terms of only two-nucleon contact terms. Pionless EFT for 
the weak NN interaction in Ref. [l| is parametrized by seven independent PV low energy 
constants (LECs) at the leading order (LO), but recently it has been shown that two terms 
are redundant and thus five terms are independent in practice {|. Thus, inasmuch as the 
number of unknown weak parameters is concerned, Danilov's idea in the past gives the same 
parametrization to the pionless EFT at leading order. 

Parity-violating vertex in the pionless theory in [l| consists of the multiplication of two 
two-nucleon fields: one in S state and the other in P state. Given a rule to transform a 
two-nucleon state to the corresponding dibaryon field, it is straightforward to obtain the 
PV Lagrangian that describes the weak A A interaction in terms of PV dibaryon-nucleon- 
nucleon (dNN) or dibaryon-dibaryon vertices. Introducing a dibaryon field for the deuteron, 
the effective range contribution ('ypd ~ 0.4) to the deuteron propagator is taken into account 
up to infinite order, and it consequently makes the convergence of the theory improved com- 
pared to the pionless EFT that does not have dibaryon fields. Since scattering lengths and 
effective ranges in the S state are unusually large, resummation of effective range contribu- 
tion to infinite order in dibaryon formalism is especially useful for the two-nucleon systems 
dominated by S state. In this work, we obtain the PV Lagrangian with dibaryon fields by 
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transforming the two-nucleon S states to the corresponding dibaryon fields, while describing 
the P states in terms of the two-nucleon fields. Weak NN interaction is described by the 
PV dNN vertices, which have unknown weak coupling constants. 

We plug the Lagrangians in the calculation of the PV polarization (P 7 ) in np — > dj at 
threshold. PV polarization has been calculated with the weak one-meson-exchange (OME) 
potentials (conventionally referred to as DDH potential [lo|]) and with various strong inter- 



action models 



The results in 11[ show strong dependence on the strong interaction 



model, and are dominated by the p- and cu-meson exchange terms in the DDH potential. In 
the EFT, p, to and heavier mesons are integrated out because their masses are very large 
scales at low-energy few-body processes, and their contributions are embedded in the NN 
contact terms. Since the PV polarization in np — > dj is dominated by the heavy mesons 
in the OME picture for the weak NN interaction, if it is considered in the EFT, only the 
contact terms are relevant and thus the pionless EFT may be one of the most favorable 
frameworks for the investigation. Result for P 7 is obtained in terms of the unknown weak 
dNN coupling constants, which have to be determined from the measurements for the rele- 
vant PV observables. 

We outline the paper as follows. In Sec. II, we present the parity-conserving and the 
parity-violating Lagrangians that contribute to the observable at leading order. In Sec. Ill, 
we obtain the PV polarization in unpolarized neutron capture by a proton at threshold, and 
discuss the result. We conclude the paper in Sec. IV. 



II. EFFECTIVE LAGRANGIAN 



Parity-conserving (PC) Lagrangian includes strong and electromagnetic (EM) interac- 
tions. PC Lagrangian with dibaryon fields can be written as 

Cpc = Cn + C s + C t + C st , (1) 

where Cn, £ s , C-t and C st represent PC interactions for nucleons, dibaryon in 1 Sq state, 
dibaryon in 3 S± state, and EM transition between 1 S'o and 3 S± states, respectively. Retaining 
the terms that are relevant to the quantity of interest in this work, we have 

C N = N* (iv ■ D + JL_ {(v Df - D 2 }) N, (2) 
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Cs = a s sl ^vD + ^- [(v ■ Df - D 2 } + A s } s a - y s {sl[N T P^ N] + h.c.} , (3) 
C t = o t t\ \iv ■ D + JL [(„ ■ Df - D 2 } + A t j ti - y t {t\[N T P^ N] + h.c.} , (4) 
C st = ^=[t\s 3 Bi + h.c], (5) 

where the projection operators for the 1 S'o and 3 Si states are defined respectively as 

^i l50) = ^<W a , (6) 
^ (35l) = ^<W*. (7) 

Velocity vector t> M satisfies v 2 — 1, and = <9 M — iV^ xt where V^ xt represents the external 

vector field. Dibaryon fields in 1 Sq and 3 Si states are denoted by s a and U, respectively, and 

— » — * 

B>i is external magnetic field given by B = V x V . a s and a t are the sign factors having a 
value —1, and A S)t are defined by the mass difference between the dibaryon and two nucleon 
states as A S;t = m s>t — 2m N . Low energy constants y s and y t are the strong dNN coupling 

2^T 



constants determined from the empirical values of effective ranges. We obtain y s = — J— 
and y t = where r is the effective range in 1 S state and p rf is the effective range 

for the deuteron. LEC L\ denotes the photon-dibaryon-dibaryon coupling constants for the 
Ml transition, and it has to be determined from experiments. 
PV Lagrangian for the two nucleon system can be written as 

£pv = ^pv ( 8 ) 
A7 

where AI denotes the isospin change in the PV vertex. PV vertex changes the orbital 
angular momentum by an odd number (e.g, S <-> P). Because A(L + S + I) has to be even, 
we have A(S + 1) = 1 for the two nucleon system. Consequently we have 

£pv = £pv ^PV (9) 

Since the total angular momentum is conserved in the NN interaction, parity mixings 
allowed by the PV interaction for the lowest orbital states are 1 S <-> 3 P , and 3 Si <-> 1 P\ 
due to £p V , and 3 Si <-> 3 Pi due to £p v . In the pionless theory interaction is described 
only by the nucleon-nucleon contact terms which have undetermined LECs. In the case 
of pionless theory with dibaryon fields, we assume that a PV dNN vertex subsumes the 
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PV NN interactions. Non-relativistic P-odd and T-even Lagrangian for the neutron-proton 
system with AI = can be written as 

£ PV = 5/2 4 ^a^r.n- (v - VW + h.c. (10) 

+— 4 iV^r ' (V - V W + h.c, (11) 



where fo°* N and /i°^ N denote the weak dNN coupling constants for the parity mixing for the 
1 S'o and 3 Si states, respectively. Spin and isospin operators o^o^Ta in Eq. (TTUT) projects two- 
nucleon system to 3 Po state. PV vertex given by Eq. ( jTUl) therefore generates 3 Po admixture 
in the 1 So state. Similarly, cr 2 r 2 in Eq. (TTTj) is the projection operator for x Pi state, and 
thus the Lagrangian mixes 1 Pi state in the 3 Si state. For the AI = 1 part, we have 3 Pi 
admixture to the 3 Si state, so the Lagrangian reads 

4v = i ^ NN 5/2 4 N T a 2 a 3 T 2 Tz % - ( V - v) iV + h.c. (12) 



Lagrangians given in Eqs. (I10llll|12p represent weak interactions between a neutron and 
a proton. Full LO interactions in the pionless theory, which include nn and pp weak interac- 
tions as well as the np one can be found in the literature [3]. By transforming a two-nucleon 
field in 5* state to the corresponding dibaryon field, one can easily get a mapping between 
pionless theories with and without the dibaryon fields. We will discuss the relation of the 
two theories in more detail in the discussion of the result. 



III. RESULT AND DISCUSSION 



In the pionless theory, expansion parameters are Q/m v or Q/A, where Q is a small 
momentum, the pion mass and A a symmetry breaking scale. Since the scattering 
lengths and effective ranges in the 1 Sq and 3 5*i states are large, we count their inverse as 
small scales, i.e. (7, l/a s , l/a t , l/r , l/pd) ~ Q, where a s u\ is the scattering length in 
1 5'o( 3 5'i) state and 7 = y/m^B with B the deuteron binding energy. Nucleon and dibaryon 
propagators are counted as 1/Q 2 and a loop integral contributes an order of Q 5 . 

Feynman diagrams at leading order are depicted in Fig. [TJ which are of the order of Q°. 
Single solid and wavy lines represent nucleon and photon fields, respectively. Double line 
with filled circle denotes the dressed dibaryon fields, which includes the infinite sum of the 
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FIG. 1: Leading order (Q ) PV diagrams for np capture. Single solid line denotes a nucleon, wavy 
line a photon, and a double line with a filled circle stands for dressed dibaryon propagator. Circle 
with a cross represents a PV dNN vertex. 

intermediate nucleon loops. Small dot at the dibaryon-nucleon-nucleon vertex denotes the 
strong dNN coupling, which is proportional to y s or y t , and the circle with a cross denotes 
the PV dNN vertex proportional to h% N , h% N and h\ NN . For the photon- nucleon coupling 
in Fig. [I] (a-c), we employ the vertex function of the convection current given by 

iT VNN (El) = ^—{1 + r 3 )~(p + p) ■ e* 7 , (13) 

where p and pi are the in-coming and out-going nucleon momentum at the photon-nucleon 
vertex, respectively, and is the polarization of out-going photons. For the PV photon- 
dibaryon- nucleon- nucleon (VdNN) vertex in Fig. [1] (d-f), we assume minimal coupling to 
the PV dNN vertex, 

V-z|(l+r 3 )y, (14) 

where V denotes the external photon field. With the minimal coupling, coupling constants 
at the VdNN vertices are the same with those at the dNN ones. Resulting amplitudes are, 
therefore, proportional to the weak dNN coupling constants h% N , or h\ NN , and thus 
we have three unknown coefficients in the result. 
PV polarization P 7 in np d'y is defined as 

where a + and cx_ are the total cross section for the photons with right and left helicity, 



respectively. P 7 was measured in 70's, and the reported value is P 7 = (1.8 ± 1.8) x 10 7 [131 ] . 



but there was no more measurement after that. At threshold, PV asymmetry in dj —>■ np 



is equal to PV polarization in np — *> dj. PV asymmetry in 



calculated with the DDH potential up to about 10 MeV above threshold 



np has been recently 
,Q. Absolute 

value of the asymmetry is maximum at threshold and it decreases very quickly as the energy 
increases. Measurement may be most feasible at threshold, and if the measurement is 
performed at threshold, it can be directly related to the PV polarization in np — *> d'y. 
Transition amplitude that includes both PC and PV contributions can be written as 



iM, 



up 



Y? d ■ (k x ?) - iZ% ■ e* 



N'P- 



es ) 



-7/ "~ c d 

Y denotes the PC amplitude, and we take the result in Ref. 161 ]. 



N. 



Y 



7 



m 



N 



1-lPd 



1 + k v )(1 - 7a s ) - 7 a s Lt 



(16) 



(17) 



where Ky (= 3.706) is the isovector anomalous magnetic moment of the nucleon, 7 (= 
y/m^B = 45.7 MeV) is the deuteron momentum, p d (= 1.764 fm) is the deuteron effective 
range, and a s (= —23.732 fm) is the neutron-proton scattering length in the 1 5'o state. We 
can reproduce the neutron-proton capture cross section at threshold, <7 ex p = 334.2 ± 0.5 mb 



16J. Z is the PV amplitude for the transition from initial So 

as 

(18) 



with La = -4.427 ±0.015 fm 

to final 3 S\ states. PV polarization P 1 is obtained in terms of PC and PV amplitudes as 

Re(YZ* 



P-, 



Y 



We obtain the PV amplitudes for the diagrams in Fig. [T] as 

,2 



dNN 



7 



P 



3m 2 N ^/m N p d \] 1 - 7p d 7 2 +p 2 ' 
1 h 0s 
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7 3 + ip 3 



3m 2 N ^/m N p d ]j 1 



1 



'YPd — s - 2 r oP +ip T+ P 
1 7 3 + ip 3 
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3m 2 Ny /m N p d 

1 h ot f 

2 m 2 Ny /m N p d y 1-jpd 

X "dNN / I 



1 - 7Pd i - l r oP 2 + ip I 2 + p 2 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 

2 m 2 N ^m N p d \/ 1 - 7p d ^ - \r Q p 2 + ip ' 

where tq (= 2.70 fm) is the neutron-proton effective range in the 1 Sq channel. Taking the 

limit p — > at threshold, we obtain the net PV amplitude Z, 

1 



7 



2m 2 NS Jm N p d 



1 



H 



Of 



l-7Pd^-- \r p 2 + ip' 
ip 



7 



7 



m 2 N y/m N p d ]j 1 - 7p d 



/f 



ot 



:7a< 



+ i/i 0s 

~ g'^dNI 



(25) 



and the PV polarization P 1 at LO reads 

| - ha s ) + ha s h° d l 



V irm N p d (1 + - 7a s ) - 7 2 a s Li 

= -(2.59^ N -l-01^ NN )xlO- 2 . (26) 

PV polarization turns out to depend on two weak coupling constants h^ n and h% N , and 
thus we cannot determine them uniquely from a single measurement of P 7 at threshold. In 
order to determine them unambiguously, we need more data for P 7 at energies other than 
threshold, or measurements of observables that are independent of P 7 . We will discuss this 
matter in more detail in the conclusion. 

Now we try to compare our result to the one obtained with a pionless theory where 



there is no dibaryon field We start from the pionless 



V Lagrangian in Ref. 



121. If we 



12j to a dibaryon field, we obtain 



transform two-nucleon fields in S state in Eq. (6) in Ref. 
the PV Lagrangian in the dEFT given by Eqs. (I10til2[) . We use the transformations from 
two-nucleon fields to a single dibaryon one given by 

for the 1 S and 3 Si states, respectively. Cq So ^ and Cq are the coefficients for the LO 
strong two-nucleon contact terms in the pionless theory. In the power divergence subtraction 
scheme, they are given as 

1 rriN ( 1 "\ 1 ttln ( 1 2 



where /x is the renormalization point. Substituting the transformations given by Eq. (f2"7j) 
into Eq. (6) in Ref. and comparing them with the PV dEFT Lagrangians in Eqs. (TTUl 



ITT), we obtain 



»£. = uJ^f K " /•) (C^o""" - 2C& ,ft ») , (29) 



2tt jv Va 

= 16^f^< ( 7 - - /x) CC*"^). (30) 

Inserting Eqs. (I29| 1301 to the result for PV amplitude in Eq. (1251) and assuming /x = m^, we 
obtain 

Z oc C^" 1 ^ - 0.56 (ci^o 3Po) " 2C^2 3P0) ) • (31) 



Isospin change is zero at the vertex denoted by C^ Sl Pl \ i.e. AI = 0, and thus assuming 
roughly C^ 1 " 1 ^ ~ C 



we obtain the ratio of the coefficient for AI = contribution 



AJ=0 > 

to that for AI = 2 one in Eq. (1311) approximately one half. PV polarization has been 
calculated in the hybrid scheme in Ref. |2j, where strong interaction is described by Argonne 
vl8 model (Av 18), weak interaction by the pionless EFT and the EM operator by Siegert 
theorem. The result in Ref. [2j is represented in terms of Danilov parameters. Substituting 
the relations of Danilov parameters and PV LECs in the pionless theory to the result, P 7 
reads 



P 7 (hybrid) = (-0.25Ci + 2.UC 3 + 4.18C 5 ) x 10" 



(32) 



where C\ and C3 correspond to AI = vertices and C5 to AI = 2 one. Assuming C3 ~ C\ 
and comparing the coefficients for AI = contribution to that of AI = 2 in Eq. (1521) . we 
obtain a ratio roughly one half, which is similar to our result. Similar value of the ratio has 
also been obtained from the calculation with DDH potential for the weak interaction and 



Av 18 for the strong one 



IV. CONCLUSION 



We have calculated the PV polarization in np — > d'-f at the threshold with a pionless EFT 
with dibaryon fields. Weak NN interactions are described with the PV dNN vertices, and 
the PV observable has been obtained in terms of the PV dNN coupling constants. Precise 
measurement of the observable will provide a constraint to determine the PV coupling 
constants unambiguously. 

EFT has been employed partially in the calculation of PV asymmetry in np — > d'y with the 



pionful theory 
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181 ]. For instance, in Ref. 
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181 ] . meson-exchange currents (MECs) 



are obtained up to an order, but the strong interaction is described by a phenomenological 
model, Av 18, and weak interaction by the DDH potential. In Refs. , on the other 
hand, weak potential is expanded up to a given order, while EM operator is accounted for 
with Siegert theorem and strong interaction described with Avl8. Calculation where EFT is 
partially employed is called hybrid calculation. Current conservation for a given PV potential 
and the corresponding MEC has been used as a crucial criterion in the calculation of the 



anapole moment of the deuteron 13, |20j. Since the current conservation can be satisfied 
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when a potential and corresponding MECs are taken into account consistently, consistent 
expansion of strong potentials, weak potentials and transition operators is an important 
requirement in the EFT. It has been pointed out that the orders of the interactions and 



transition operators in the hybrid calculations are in serious disagreement [12]. In our 
calculation with pionless dEFT, the order of a diagram is obtained by counting the strong, 
weak and EM vertices altogether, and we truncate the expansion at a given order. Therefore, 
our calculation satisfies the consistency requirement mentioned above. On the other hand, 
results from the conventional calculation, where strong interaction is accounted with modern 
potential, weak interaction with DDH potential, and EM operator with Siegert theorem, 
have provided benchmarks to both experiment and theory, but the physical criteria such 
as current conservation have seldom been checked carefully. It is important to understand 
the uncertainty due to the order mismatch in the conventional and hybrid calculations, and 
investigation along this direction with either pionful or pionless theories is an important 
future work. 

There are five weak LECs in the pionless theory and therefore we need at least five data for 
the PV observables. P 7 may be one of them. Recently PV longitudinal asymmetries in pp, np 
and nn scattering have been calculated with a pionless EFT [12]. Longitudinal asymmetry 
in pp depends on three PV coupling constants Ca/=o 12°\ an d thus the measurement at 13.6 
MeV provides a relation for them. PV asymmetry in np — > d'-f and deuteron anapole moment 
have been calculated with the pionless dEFT 2l|], and the results turn out to be dominated 
by h\ NN . Measurement of the PV asymmetry at SNS is expected to play an important role 
in determining the value of weak LEC hl NN (or the weak pion-nucleon coupling constant h\). 
Turning to the possibilities in the three body system, one can find a recent calculation of the 
weak effect in the spin rotation in nd scattering 22j. The authors employed DDH potential 
for the weak interaction, and obtained a result dominated by h\. Though the asymmetry in 
np and spin rotation in nd are observables independent to each other, they are exclusively 
dependent on h\ (or equivalently h\ NN ) , so the measurements of the observables will provide 
a check for the consistency of h\. In order to determine the remaining weak LECs in the 
pionless EFT, we need calculations and measurements for as many observables as possible. 

Among many possible PV observables in the few-body systems, an observable that draws 
our interest is the PV asymmetry in nd — > £7 at threshold. It has been measured at ILL 
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231 ] . and the reported result reads 

A\ = (4.2 ±3.8) x 10" 



Theoretical calculation of in Ref. [24| adopted DDH potential for the weak interaction, 
and examined the dependence on the strong interaction models such as de Tourreil-Sprung 
(TS) and Reid soft core (RSC). The results are interesting in some aspects. First, dependence 
on the strong interaction model is non- negligible; TS model gives a result A* (TS) = 0.81 x 
10~ 6 while RSC gives v4^(RSC) = 0.61 x 10~ 6 . Second, isoscalar, isovector and isotensor 
PV interactions in the DDH potential give similar contributions to A*, e.g. 0.40, 0.45, and 
—0.04, respectively, with the TS model. This means that contributions from tt, p and uj 
exchanges in the PV potential are similar to each other. Dependence on the strong model 
and the non-negligible contribution from the heavy mesons are the features common with 
the PV polarization in np — > d^y. In this problem again, therefore, pionless EFT will serve us 
with a most natural and systematic way to parametrize the parity-mixing in the few-body 
system due to weak nuclear force. 

Parity violation in the three- and few-body systems can show us the effects that are not 
accessible in the two-body systems. For instance, strong 3N force can give non-negligible 
correction to the one- and two-body contributions to the PV observables. There has been 
no consideration on the weak three-body force, but we have recently obtained non-zero 



25j. Two- and three-body PV 



component of weak 3N force in a preliminary calculation 
meson-exchange currents are also important issues. We expect that the EFT will play a 
crucial role in extending our understanding of the nuclear weak force in few-body systems. 
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